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A complete study of the one-loop induced decay Z — > uwy is presented within the framework of 
the Standard Model. The advantages of using a nonlinear gauge are stressed. We have found that the 
main contributions come from the electric dipole and the magnetic dipole transitions of the Z gauge 
boson and the neutrino, respectively. We obtain a branching ratio B(Z — > vwy) — 7.16 x 10 -10 , 
which is about four orders of magnitude smaller than the bound recently obtained by the L3 Collab- 
oration and thus it leaves open a window to search for new physics effects in single-photon decays 
of the Z boson. 
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I. INTRODUCTION 



In the absence of any clear deviation from the Standard Model (SM) predictions, considerable work has been done 
recently to search for processes that are forbidden or highly suppressed and may become a good playground to test 
any new physics lying beyond the SM. The radiative decay Z — ► P1/7 is a good example of this type of processes 
(T) ' because it arises at the one-loop level in the SM. The L3 Collaboration has searched for energetic photon events 
coming from single-photon Z decays. It was obtained a bound of one part in a million for the respective branching 
ratio Single photon events in Z decays have been associated to gauge boson and leptonic substructure ^,|), 
supersymmetric particles Q, a possible neutrino magnetic moment || and the production of a Higgs boson with 
' its associated CP-odd majoron 0. The decay Z — > vv^j is expected also to be sensitive to the possible anomalous 
electromagnetic properties of the Z boson and the neutrinos. The latter point has become of increased interest after 
O ' the recent evidence of oscillations in atmospheric neutrinos . In a recent study of the decay Z — > vv^( within the 
effective Lagrangian approach |^,^J , the L3 bound was used to get direct constraints on dimension-six and dimension- 
' eight operators associated with the anomalous electromagnetic properties of the neutrino. 



Though there is a great incentive in studying single-photon Z decays, to our knowledge its decay width has not 
been computed in the SM. We would like to stress that in order to disentangle new contributions to this process, 
it is important to count with a complete analysis of the SM contributions. Our general aim is precisely to present 
explicit expressions of a calculation of the SM decay Z — > vv^j at lowest order in perturbation theory. As we will 
see below, the amplitude for this process can be written in terms of five independent U e (l) gauge structures. Some 
of them can be identified in terms of the so-called off-shell electromagnetic properties of the Z gauge boson and the 
neutrino. It is well known that the Z boson possesses both an anapole moment (AM) and an electric dipole transition 
(EDT) defined through the ZZ7* and ZZ*^ couplings, respectively (where * stands for off-shell particles) A SM 
neutrino can have only off-shell electromagnetic properties characterized by the vvy* and vv*^ couplings. The vv^j* 
coupling corresponds to an AM of the neutrino, which has been extensively studied in the literature [jl(J[ll). It will 
be seen below that the w*j coupling can be identified with a magnetic dipole transition (MDT) of the neutrino. In 
this paper we have found that the main contributions to the decay Z — > z/z/7 come precisely from the MDT and EDT 
properties of the neutrino and the Z boson, respectively. In particular, the MDT of the neutrino plays an important 
role in this process. Therefore, we will discuss with some extent its SU(2) structure. For this purpose we will use 
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a nonlinear gauge which allows us to separate in a transparent way the SU(2) and U e (l) gauge structures and, as a 
consequence, it will simplify considerably the calculation. In particular, we have found that there is no SU(2) gauge 
link between the coupling and the box diagrams contributing to the decay. This result should be contrasted with 
the situation found in the calculation of the SM scattering amplitude of neutrinos with charged leptons and nucleons. 
In this case, the SU(2) gauge structure of the 1/1/7* coupling is directly linked to the box diagrams contributing to 
the scattering process fLof . 

In the nonlinear gauge the decay Z — > Pwy receives contributions from reducible triangle diagrams characterized 
by the ZZ*~f and vv*"j couplings, as well as from box diagrams containing W boson and lepton combined effects. It 
is well known that in many four-body processes the box diagrams contributions can be neglected compared with the 
triangle ones. However, this approximation fails when the triangle and box diagrams possess common U c (l) gauge 
structures which induce form factors of the same order of magnitude. In particular, the tensor structures arising from 
the box diagrams will have a Lorentz covariant decomposition inducing the MDT of the neutrino. Therefore, it is 
not possible to neglect the box diagram contributions in the decay Z — > Dvy. Our numerical result confirms that this 
decay is a SM highly suppressed process, about four orders of magnitude below the L3 bound 0. Thus it leaves open 
a window to search for new physics effects in single-photon events, mainly in the neutrino sector. 

We have organized our presentation in the following way. In Sec. [n] we present the analytic expressions for the 
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Z — > vvy decay width. In Sec. Ill we give details of the numerical analysis, and the conclusions are presented in Sec. 



II. THE NONLINEAR GAUGE AND THE DECAY Z -> P1/7 

The decay Z — > Pz/7 is a SM one-loop prediction which receives contributions from all charged particles. Within 
the conventional linear gauge, we can classify in five sets the Feynman diagrams contributing to this decay at the 
one-loop level. The first three sets are characterized by reducible diagrams including the Z* — 7 mixing term (Fig. 
|l|) and the ZZ*^ and vv*^ couplings (Fig. || and Fig. ||, respectively), while the other two sets include box diagrams 
where the Z boson couples to charged leptons (Fig. ^) and to W bosons (Fig. ||). For the purpose of our calculation, 
it is suitable using a nonlinear gauge in which there is no SU(2) gauge (nor U c (l) gauge) link between the vv*^ cou- 
pling and the box diagrams. This task is a difficult one to fulfill in a linear gauge, where SU(2) gauge independence 
and U e (l) gauge invariance are only accomplished after adding up almost all the Feynman diagrams involved in the 
calculation. Furthermore, in this type of gauges a renormalization of the Z* — 7 mixing term has to be done In 
contrast, it can be shown that in the nonlinear gauge the amplitude for this mixing term vanishes for any value of the 
SU(2) gauge parameter £. We will show also that the vv*') coupling is by itself SU(2) gauge independent and U e (l) 
gauge invariant. 

A nonlinear i?j-gauge condition was first introduced by Fujikawa [fl3|| . The ordinary derivative involved in the 
gauge-fixing term for the W boson is replaced by the respective U e (l) covariant derivative; as a consequence, there 
are no VK ± G T 7 unphysical vertices, with G ± the respective would-be Goldstone boson. The W boson and the charged 
ghosts sectors are separately U c (l) gauge invariant, leading to naive Ward identities. This procedure was later ex- 
tended to remove both W^G^j and W^G^Z vertices [ [l4|Jl5| . It is possible to proceed further with this scheme and 
remove more unphysical vertices, such as H°W ± G T j(Z) and G°W G^j(Z), with G° and H° the neutral would-be 
Goldstone and physical Higgs bosons, respectively For our present purposes, we will use a gauge-fixing func- 

tional which is nonlinear in the vector sector but linear in the scalar sector. 

The functionals which define the respective W, Z and A propagators in the nonlinear gauge can be written as 

/+ = DpW+i* - itm w G + , 
f z = d li Z»-£m z G , 

f A = M", (i) 

where s w — sin9 Wl c w — cos8 w and 8 W is the weak mixing angle. We have also defined the following operator 

D [j, = Op ig'Bp, (2) 
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Z^ = c u>W p — s w B^, — s w W^ + c w B^, and W p and B p are the gauge fields of the SU(2) and Uy(l) groups, 
respectively. We can see that the above operator contains the electromagnetic covariant derivative, so that the / ± 
functionals transform covariantly under the U (l) group. The corresponding gauge fixing Lagrangian is given by 

Car = -\ f + r ^ (f Z ) 2 ^ (f A ) 2 , (3) 

which removes the W^G^'y and W^G^Z vertices from the Higgs kinetic energy term {D ^(fy [D^ (f) . The correspond- 
ing Faddeev-Popov Lagrangian is given by 

C FP = -c-[D fl D' i + £m w (m w + H° + iG°)]c + - igc w c- D p {W + ^c z ) - iec~ D^W^Cy) 

2 

_ i^EW + ^c-(d"c z ) + ieW+c-{d^ Cl ) - igc w W + ^(d^c z )c- - ieW+^id^c' 
— gc 2w mz£,{G + c~ cz + G + czc~) — eniw£,{G + c~ c 7 + G + c 7 c~) + H.C. 

-c z [a+£mz(mz + H )}cz-c 7 ac 7 , (4) 

where c ± {c ± ) 1 cz(cz), and c 7 (c 7 ) are the pairs of ghosts associated with the W ± 1 Z, and A gauge bosons, respectively, 
and c 2w — <? w — s^. The charged ghosts satisfy (c + )^ = cT and {c + Y = c~ . We have introduced the following operator 

D p = 3^ - igWl, (5) 

which also contains the electromagnetic covariant derivative []. Therefore, just as it was the case with the f~f + term 
in Cgf, the Cfp Lagrangian is invariant under the U c (l) group. Since both the Yang-Mills and the kinetic energy 
Higgs sectors can be expressed in terms of the D p and D p operators, all the charged sector of the SM satisfies naive 
Ward identities. 

To calculate the amplitude for the decay Z — > vvj in the nonlinear gauge, the Feynman rules for the 
A x (ki)W+(k 2 )W-(k 3 ), Z x {ki)W+{k 2 )W-{k i ), and A p Z x W+Wp vertices are required 

T \™ W = _ie p3 -ki + y) P .9Ar, + {ki -k 2 - y )v9\ P + (h - k 3 ) x g PV ], 



^x PV = -~igc w [(k 3 - ki - -f-r) P g\n + (hi - fc 2 + '-f-p) v g\ P + (k 2 - k 3 ) x g P r,}, 

C l« ? C w S 

1 S 2 

r P \$ W = ~iegc w [2g af3 g pX - (1 + --f){g a xgp P + gp\g ap )}, (6) 

S c w 

where all the momenta are incoming. We have taken the lepton masses to be zero except where large logs (collinear 
singularities) arise. This approximation and the use of the nonlinear gauge reduce considerably the number of dia- 
grams involved in the calculation. We have calculated the Feynman amplitudes using the computer program FeynCalc 

In the linear gauge the Z* — 7 mixing term receives contributions from W gauge bosons, charged ghosts and 
W ± G T combined effects ^. In contrast, in the nonlinear gauge the H /± G =F combined effects are absent. From (Q), 
(0) and (|4|) we can see that the f + f~ term in the Cgf and the Cfp Lagrangians also presents explicit U (l) gauge 
invariance. Consequently, in this gauge we expect a zero contribution arising from all charged particles. This fact 



was shown before in the Feynman-t'Hooft version of the nonlinear gauge JL7 18 1. We will show that this result is 
valid for any value of the gauge parameter £. Notice that the charged part of the Cfp Lagrangian resembles scalar 
electrodynamics (see (^)), so that we do not expect anything new coming from these particles in the general nonlinear 



Notice that if we use in (hh _D M instead of _D M , as it was pointed out in Jl|, the W ± G^ 1 vertex is removed but the W ± G^ Z 
vertex only becomes modified. 

2 Since the Lagrangian which defines the G boson is not modified in its interactions with the photon by the gauge fixing 
procedure used for the SU(2) group, and since it resembles scalar electrodynamics, the contribution of this field to the Z* — 7 
mixing term vanishes in any gauge. 



3 



i?£-gauge. In contrast, both the I^-propagator and the W ± W T, -f(Z) vertices have a nontrivial dependence on the £ 
parameter. The corresponding amplitude which arises from the diagrams shown in Figs. [fj(a) and [j](b) is given by 

n 5 7 = i -^^[u a + u b + |(K + H)]fa* (7) 

with U a (Ub) and V a (Vb) characterizing the loops in Figs. |l|(a) (0(b)). After a straightforward calculation, we have 
found that there is not Z* — 7 mixing term in the general nonlinear i?£-gauge 

U a = -U b = ei r o + 6B (0, £m 2 w , tm 2 w )] - 3[1 + 6B (0, m 2 w , m 2 w )}, 

V a = -V b = -^[3 + 2B (0, Zm 2 w , £m 2 w )} - ~[5 + 6S (0, m 2 w , m 2 w )}, (8) 

where Bq(0, ^m^, £to^/) and Bq(0, m^, m^) are Passarino-Veltman two-point scalar functions |2C|| . 
In this manner, the width for the decay Z — > Pz/7 can be expressed as 

r l r l-x 

m z [ - f ,.1,-^12 



r(Z -> ^7) = y o J dy|A4| J , (9) 

where is the invariant amplitude 

M = uMMf, +M a J + Mf ox4 + Mf ox5 ]v R (j> 2 )e a (k 2 ,\ 2 )e* (ki, Ax), (10) 

and the M. a z », , M^ oxi , M^f ox5 amplitudes correspond to the respective diagrams shown in Figs. ||~||. The 

explicit form of these amplitudes will be given below. Moreover, pi(p 2 ) are the neutrino (antineutrino) momenta, 
ki(k 2 ) are the "f(Z) momenta, ep(e a ) are the "f(Z) polarization vectors, and tf£=|(l — 75)11. It is useful to define the 
following dimcnsionless quantities 

x = o — , y = 7j — . (11) 

and the combinations 5 = 1 — x — y, 6+ = x + y, 8- = x — y, S x = 1 — x, 5 y = 1 — y. 



A. ZZ**y coupling contribution 

The ZZ*"f coupling contributes to the decay Z — > vv^i through the reducible diagram shown in Fig. |^. The 
contributions to this coupling are due to the fermion triangle diagram responsible for the ABJ anomaly, which of 
course vanishes with the SM fermion multiplet assignment. The corresponding amplitude must be U (l) gauge 
invariant by itself. This amplitude was obtained by Barroso et al. some time ago || . We have reproduced this result 
to compute the corresponding contribution to the decay Z — > vvy. Since one Z boson and the photon are on-shell, 
we take k^ep = k^t a = 0, k\ = 0, k\ — m|. In addition, using Bose statistics and the Ward identity satisfied by the 
electromagnetic current, we obtain for this amplitude 

= jS* zZvWaQi { -T^T^jCoirriU 2 ) 

%^|) (jB ° (m|) ~ Mq2)) + \^ XPk2 ^ 

+ \i^^{ m U 2 ) ~ -f^r(B (m%) - B (q 2 )) + ±]e^*u}, (12) 

Zi y TJl g y Til g Zi 

where C (m 2 z ,q 2 ) = C (0,m 2 z ,q 2 ,m 2 ,m 2 ,m 2 ), B (m 2 z ) = B (m 2 z , m 2 , mj) , and B (q 2 ) = B {q 2 ,m 2 ,m 2 ) are 
Passarino-Veltman three- and two-point scalar functions written in the notation of . Q f is the electric charge in 
units of e, q 2 — (pi + P2) 2 = Tn z 5, g v = £3 — 2Q-^s 2 v , and g^ A = t^, with t$ the third component of weak isospin. 

Now, it is straightforward to obtain the respective contribution to the decay Z — » z/1/7 
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Mf, = J^ lx Y,9v9 f A Q f F(m 2 z , q 2 ,rn})e^k lp , (13) 

yS w C w ) 

where 



1 fP'Tfi 

F{m%,q\m}) = [(g 2 + m%)m}C (m%, q 2 ) «^(5 (m|) - B ( q 2 ))}. (14) 

[q m z) 9 m z 

In addition, we have removed from ( |i"4| ) the mass independent term which vanish when we perform the summation 
over fcrmion families. 



As it was shown in ||, and more recently in pl| , the Z boson possesses only off-shell electromagnetic properties 
characterized by an electric dipole transition (EDT) and an anapole moment (AM), which correspond to the ZZ*<y and 
ZZj* couplings, respectively. Since the on-shell ZZj vertex is zero, the Z boson does not possess on-shell (physical) 
electromagnetic properties. So the contribution of the ZZ*j coupling to the decay Z — > i>i/j corresponds to an EDT, 
which is a transfer-momentum dependent quantity. This EDT electromagnetic property of the Z boson depends only 
on the number of families. Indeed, as we will discuss in the following section, the contributions arising from the first 
two families, which are nearly mass degenerated, are very small compared with the corresponding contribution of the 
third family, where the masses involved have an important gap. 



B. coupling contribution 



The decay Z — > vv^j receives contributions coming from the vi>*j coupling through the diagrams shown in Figs. 
||(a) and [|(b), where the photon couples to leptons and to W bosons, respectively. We will present the calculation for 
the vv*^f coupling in the general nonlinear i?£-gauge and then we will use it to give the corresponding amplitude for 
the decay Z — > vv^. The triangle ||(a) depends on the gauge parameter £ only through the W-propagator. On the 
other hand, the triangle ||(b) has a more complicated £ dependence through both the W-propagator and the W^W^j 
vertex. Taking into account that a neutrino and the photon are on-shell, the amplitude for the vi>*j coupling can be 
written as 

r 7* 7 = iCl^Fa + F b ) +px4x{G a + G b )], (15) 

M I67H p z 2 

where p 2 = (fci + pi) 2 = xrivi is the momentum transferred by the Z boson. F a (Fb) and G a (Gb) are the amplitudes 
corresponding to diagrams ||(a) (|](b)) 

K = ^ + 5(1 - + |d - - ^)Bo(p 2 Am 2 w ) 

- ^B {0,Zm 2 w ,Zm 2 w ) + h^-Z)B (p 2 ,0^m 2 w ) 

Z Z TTl^r 

3 m 2 

-(-- ^f)S (0, m%r,m%r), (16) 
z p 



G„ = l 



2(1 



>w 



)[B Q (p 2 ,0,m 2 w )-B (0,in 



(17) 



F h 



2 { m 2 



;(1-0-tt(1 



^)(1 



nv 



)B (p 2 , 0, m 2 w ) + -£Bo(0, £,m 2 w , ^m 2 w ) 



w 



-e)Sob 2 ,0,^) + (|-4 
2 p z 



)Bq(0, mlv, m 2 w ) + 2{m 2 w - p 2 )C (3), 



(18) 
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G b = -3+ — f- + 2(1 - -f )[B (;p 2 , 0, m 2 w ) - B o (0, m 2 w , m 2 w )\ - 2(m 2 w - p 2 )C (3), (19) 

where the three-point scalar function is defined in the Appendix. From (|l5|)-(|l9|) it is evident that the amplitude for 
the vv*^f coupling is a SU(2) gauge independent and U e (l) gauge invariant quantity that can be written as 

T up = i A(p_ k , a ^ (20) 

where = | (7^79 - 7/37/J, and 

A(p 2 ) = ^[l + (m 2 w -p 2 )C (3)}. (21) 

It is well known that in the SM the neutrino have no static electromagnetic properties. It is easy to see that our 
result ( |2l| ) is consistent with this fact by noting that the A(p 2 ) form factor vanishes when the neutrino is on-shell. 
Using the above results, it is straightforward to calculate the corresponding amplitude for the decay Z — > P1/7. Taking 
into account all diagrams shown in Fig. [|, we have 

= o ^ 2 [Ai(x)Tf+A 2 (y)Tf}, (22) 

where the form factors A\ (x) and A2 (y) are transfer-momentum dependent quantities given by 

A 1 (x) = -[l+m%(cl-x)C (3)], A 2 (y) = A 1 {x — » y). (23) 

x 

and the U c (l) gauge structures are given by 

Tf = lA^V, Tf = ik ltll a a^. (24) 

The M a f, amplitude can be interpreted in terms of the off-shell electromagnetic properties of the neutrino. In fact, 
we can identify in the gauge structure of this amplitude the familiar form of a fermionic magnetic dipole transition. 
Since the corresponding form factor is an on shell vanishing quantity, we can relate it with a MDT of the neutrino. 
We would like also to stress an aspect related with the off-shell electromagnetic properties of the neutrino, namely, 
the MDT an AM. We have found that the MDT is an SU(2) gauge independent and U c (l) gauge invariant quantity. 
This result must be contrasted with the one obtained for the AM arising from the vvf coupling which depends on 
the gauge parameter £ • 



C. Box diagrams contribution 



In this section we present the amplitudes for the box diagrams contributing to the decay Z — > vv~f. These diagrams, 
shown in Figs. |J and|^, induce terms with the same Lorentz structure of the neutrino MDT. This implies that we can 
not neglect these contributions because they are of the same order of magnitude than the one induced by the vv*^( 
coupling. We have found that there is not Z* — 7 mixing term in the general nonlinear i?^-gauge. Also, we have found 
that the vv*^f coupling is SU(2) gauge independent and U e (l) gauge invariant. This implies that the contributions 
arising from the box diagrams must respect these properties by themselves. Even more, since the two sets of box 
diagrams arise from different Z boson couplings, we expect that they are separately SU(2) gauge independent and 
U e (l) gauge invariant. We have found that in the Feynman-t'Hooft version of the nonlinear i?^-gauge the respective 
amplitude is U e (l) invariant and finite for each set of diagrams 

M ^ = ^^t Al ^f (25) 

and 



G 



2 7 

a c„ 



^i 5 = -^E^(^ y )rf, 

s ™ m z i= i 

where the rf and r f U (l) gauge structures are given in (|24|) and 

7 Q (fci -pipf - fci -p 2 Pi), 

7/*P?(*i ■ PiS"" " rf = rf (p x <- ft), 

while the A\ and functions are given by 

4(*>tf) = ^^{5-150(5) - fl (6)] - s<y+[B (6) - S (l)] - I/[Bo(6) - B (2)]} 
2xyo+ 

{5-5+[x(x + 28 y )y 2 + 2xy(l + x)c 2 w + S-8 + ct,)C (3) 



r af3 _ 
1 3 — 


1 




m| 


r a/3 
1 4 — 


1 




m| 


pQ/3 

1 5 — 






m| 



4<5x 2 y 2 

+ y<5[x(5a; + 2S y )y 2 + 2xy(2S+ - S x )c 2 w + 5-6+<*]Co(4) 
+ a:[arf(a;J - 2y 2 ) - 2x(x - 5_ ( 5 + )c 2 t! - 8-5+0% ]C (7) 
+ 2/[<5 2 2/ 2 - 2 2/ (x<S :c + y^)c 2 u - a_*+c* ]C (8) 
+ z[y<5 x (5 2 + 2(4 + 2y)) + 2(<5 2 + xyS - yS y )c 2 JC Q (9) 

+ yS y {yS y (4x - 5 y ) + 2(x - y S y )c 2 w }C (lO) + 8 2 {SSS + + 2xy 2 + S-6 + c 2 JC {ll)} 

4 

+ IfeflF ^ + 2S v) x2 y 3 ~ xy 2 (x(3$ + 10) - ZySyK + xy{5 2 + 2x(l + 3y) + 45 y )ci 

- SS+c^Doil) + (xS + 5 + c 2 w )[x5(xS - 2y 2 ) + 2x(x - 5-8+)<* + 8-8+<*]Do(ti) 

- (yS + S + cl)[y 2 S 2 - 2y(6 2 + - y)c 2 w + 6S +t *]D {6)}, 

A l 2 (x,y) = -A\(x^y), 
A l 3 (x,y) = ^{[B (6) - B (5)] + 2xy[B (l) - B„(6)]} 

:{x[xy(3xy + 2S 2 ) + 4xyc 2 w + 8 2 c 2 w ]C (3) 



2x 2 y 2 5 2 

+ x[x(x + 2y)8 2 - 2x(x 2 + y(y - 28 + ))c 2 w + 8 2 + ci]C Q (7) 

- x[x 2 {3y 2 8 x + y 3 ) + 2x(y 2 + xS 2 )]C (9) + x 2 8 2 (5 + 2c 2 )C (11)} 

{[]-xy{8 2 - 3xy) + (2x 2 (x 2) + xy + l -8{2 + 3xy))c 2 w + h(S 6)4]D (4) 



2x 2 y 2 5 2 u 2 ay a ' y v a 2 v *" w 2 

- 2 + c5 2 c 4 



(xS + 8+c 2 J[x8 2 (x + 2y) + 2x(S x (x + 2y) - y 2 )c 2 w + 8 2 + 4]X> (5)} + (x ~ y), 



2 

+ ^2 - 2 ^ 2 ) + y4]Co(3) + y{x8 2 x + yc 2 y )C (4) 

+ xyc 4 w C (7) + y 2 (8 + c w ) 2 C (S) - x8lc 2 w C (9) 



+ [-y^5 2 + tf„(x - 25_<5 + + 2y^)c 2 u + -=-(y(2 + y) - l)c 2 ]C (10) 

+ y( S 2 (5 + 2 c 2)C (ll)} 
4 

+ ^2 - + ^2/(2 + a:) - y(3 - y))c 2 w + yS + ci}D (4) 

+ y(x5 + 8+cl)[x(x + 2y)8 2 + 2x(S x (x + 2y) - y 2 )c 2 w + 8 2 + c 4 w ]D (5) + y(y8 + 8 + c 2 J(8 + c 2 u ) 2 J D (6)}, 
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A l 6 (x,y) = J- + ^ {fa/Ms - 25„) - ^ + x 2 d+}[B Q (6) - B (5)] + x 2 S 2 + [B (6) - B (l)} 
+ yS 2 + (x-26 v )[B (2)-B (6)]} 

„2„,2 



{[x(-Sx 2 y 2 - 2xy(x(2x + 5 X )) - 5)c z w - (x6_6+ + 2y 2 S)c%]C {3) 



2x 3 y 2 S 2 

- [3x 3 y 2 + 2xy(55 x + 2x 2 )c 2 w + {xS_S+ + 2y 2 8)c 4 w ]C (A) 

+ x[-x 3 6 2 + 2x(x{6-6+ - x) + 2y 2 S)c 2 w - (xSS+ + 2y 2 5)c 4 w ]C {7) 
+ y[y 2 5 2 (x - 25) + 2y(x 2 S + - y(xS - 25 2 ))c 2 w - {x8.5+ + 2y 3 S)c 4 w ]C (S) 
+ x[x 2 6 x (3y 3 + S 2 ) + 2(x 2 (yS x + S 2 ) + y5(l - 2x))c 2 w ]C (9) 
+ y6y[y(x(3x 2 - 5 2 ) + 2S 3 ) + 2(x(x 2 - S 2 ) + 2yS 2 )c 2 w }C (W) 

- 6 2 {6(x5_5 + + 2y 2 S) + 2(xyS + + 2y 2 6)c 2 w ]C a (ll)} 

4 

+ 2a ,r^ 2 {[ 3 ^V + x 2 y 2 (-25(y + 35) + x(5y + 66) + 7x 2 )c 2 w 

+ xy(6yS 2 + x(y(y - A5) - 2S 2 ) + 2x 2 (3y + 5) + 5x 3 )c 4 w + (x5S + + 2y 2 5)c 6 jD (4) 

+ {(xS + S + c 2 w )(x 3 S 2 + 2x(x(x8 x + iy 2 ) - 2y 2 5 y ))c 2 w + (2yS y + x(x 2 - 2,y 2 ))c 4 w ]D a {h) 

+ [y 2 5 2 (2S y - 3x) + 2y(yS y (2S y - 5x) + x 3 )c 2 w + ((2S y - y)y 2 - 3xy 2 )c 4 w ]D (6)}, (32) 

A l 6 (x,y) = -A{(x~y), (33) 

A\{x,y) = -A l 5 (x~y), (34) 

AY{x,y) = ^-L-{5_[B (4) - Bb(3)] + S + [x(B (l) - B (4)) - y(B Q (4) - B (2))]} 

+ ^2- s ti 6S +(- x2s + y 2 ( 7x + V 5 )) - 2(56-(x 2 + y 2 ) + 4xyS 2 + )c 2 w - 26_6 2 + c 4 w }C (l) 

- 5 2 [x 2 (x - S y ) + y 2 (3x + 5 y ) + 25-5+c 2 JC (2) 

+ x[xS(x5 + 2y 2 ) - 2x{5- + 5 2 + )c 2 w + 8-5+c 4 w ]C {Z) 

- y[y 2 S(5x + 5 y ) + 2y{5-{\ + 5_) - 4z 2 )c 2 + 5-5+c A w ]C (±) 

- x[l - 3y 2 - 6+(2 - 6+) - 2{x + y(xy - S y ) -x 2 (2- 6+))c 2 w }C (5) 
+ yS y [yS y (2x + S y ) + 2(x - yS y )c 2 w ]C a (6) 

+ [-x(x + 2S y ) + 2xy(x + 2y)c 2 w + 6-5 + c 4 w }{xC (7) + yC a (8)}} 

4 

+ -^-^{[x5(3x(y - x8 x ) + 2y(2x 2 + y 2 ) - IZxy 2 ) + x 2 S(x(l + x 2 + y 2 ) 

+ 2S y (xy - 5_5 + ))c 2 w + x(S+(y(2x + 7y - 3) - 3xS x ) + 8xy)ci + 5-6%<&]D {l) 

+ [y 3 S(x(4S y -x) + S 2 y ) + ? y 2 5(35_(l + S-) - 4xy)c 2 w 

+ y(x(x(x - 7y) - 2y) - 3<5_(x + yS v ))c 4 w + 5_<^cg,]Do(2) 

+ (xy - 5 + c 2 J[xy 2 (x + 25 y ) - 2xy(x + 2y)c 2 w - 5-6 + c 4 jD Q (3)}, (35) 
A$ r (x,y) = -AY(x~y), (36) 

Af(x,y) = ^- 6 {B (3) - B (4) + 2x[B (4) - B (l)]} 

{-x 2 5 2 {5 - 2c 2 u )C (2) + x[x6(x - 5+(x + 2y)) + 2x(x5 - 2y)c 2 w - 6 2 + }C (3) 



2x 2 y 2 5 2 

+ S + [x 2 5(5 x - 5y) - 2x(xS - 2y)c 2 w + S + c 4 JC (l) 
+ x[xS x (6 x (5 -y)- 2y 2 ) - 2x(x + <5 2 )c 2 u ]C (5) 
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+ x[(xy(2(xS x + yS v ) - xy) - 2xy(l + S)c 2 w - 5 2 + c 4 w ]C {7)} 

Tft^ \ 

+ 2x 2 y 2 5 2 ^ XV ~ S+c ™^ xy (~2 Xy + 2xSx ^ ~ Xy<yl + ~ x5 + c Vl D o{ 2 >) 
+ [x 2 S 2 (xS x - y(3 - 2y)) + x5(5xyS y + y(x 2 + 2y 2 ) + 3x 2 5+)c 2 w 
+ x5+(3x6 x - 2xy + y(y - b))c 4 w - *1<&]D (1)} 

+ {x «-> y), (37) 
Af(x, y) = - J- + -^L^{2(2:r + y)[B Q (4) - B (3) + y(B (2) - B (4j)} + x 2 (B a {2) - B (3))} 

iff) 2 

+ ^-^{S y [2yS + + 2S(x 2 (S_ - 1) - y 2 )c 2 w + yS + S 2 c 4 w ]C (l) + yS y S 2 (S - 2c 2 JC (2) 

+ x8 x c 2 w [x5{2y - 8 X ) + yc 2 JC (3) + y6 x [y8 2 + S(xS x + 2y)c 2 w + yc 4 w ]C {4) 
+ x8 y 8 x c 2 w C (5) - [y88 2 y - (8 2 y (28S+ + 2y8 x - x(l + x 2 )) + 2x 3 y 2 )c 2 w ]C (6) 

- y8 y (x8 y - c 2 Jc 2 w [xC (7) + 7C (8)]} 

+ ^2^^ + xS ( xS * v( 2 x )) c l + ^+4]A,(i) 

- [S 2 y S 3 - yS 2 (x(x - 2) + 3y)c 2 w + 8(x5 x 5_ + 3y 2 )c 4 w - y8+c 6 w }D (2) 

+ y(xS x - c 2 w ){xy - 6 + c 2 w )c 2 w D (3)}, (38) 

A ^ {x ' y) = x^ + x^55\ {[5+iy{ ' iX 25y) ~ x2) 2 *o(3) - B (4)] 
+ 5\ [y(x - 2<y„)(B (4) - B (2)) + x 2 (B (l) - B (4))}} 

2 

+ ^^-^ {5+5 2 [S(x(x 2 + by 2 ) + 2y 2 5) - 2(x(x 2 + y 2 ) + 2y 2 S)c 2 w 
- 2S 2 + (25y 2 + xS-S+)4,]Co(l) 

+ [S 2 + 5 2 (x(x 2 + by 2 ) + 2y 2 5) - 2(x6 + 5_ + 2y 2 5)c 4 w + 2S + (SS+(x 2 + 3y 2 ) 
+ 2y(6 2 (x 2 +y 2 )-2x 2 6 + ))c 4 w }C a (2) 

+ x5+[x 3 5 2 + 2x(-xSS- + 2xyS+ - yS)c 2 w + (x5_5 + + 2y 2 S)c 4 w ]C (3) 

+ yS+[y 2 S 2 (5x + 25) + 2y(2x 2 S+ + 5{x - 25 x ))c 2 w + {xSS+ + 2y 2 5)c 4 w ]C a {A) 

+ xS+[x 2 5 x (3y 2 - S 2 ) + 2(x 2 (yS x + S 2 ) - yS(2x - l))c 2 w }C (5) 

+ yS + S y [y(x(3x 3 - 5<5 2 ) - 25 3 ) + 2{x 3 - 5S 2 )c 2 w }C (6) 

+ [-3x a y 2 + 2xy(2y(x -5) + x 2 )c 2 w + (x5_S + + 2y 2 6)c 4 w ][xC (7) + yC (8)}} 
M 4 

+ 2x3 ^ s2 {[-x 4 5 3 + x 2 5(5(3xS- + 2y(l - 2y)) - 4xy5_)c 2 w 

+ x{x5 2 + {4y - 35) + 35S+ - xS + 2yS x )c 4 u + 6+(xS-5+ + 2y 2 S)c 6 w ]D (l) 
+ [-y 3 S 3 (3x + 2yS) + y 2 5(4:X 2 5-5(x(5x - 3y + 4) - 6y))c 2 w + (~3xy5(2x + y) 
+ 2S 2 (3y 2 - xS-) + x 2 (A(x 2 + y 2 ) + x(7y - 5 x )))c 4 w ]D (2) 
+ (xy - c 2 w S + )[(3x 3 y 2 - 2xy(2y(2x - S y ) + x 2 ))c 2 w 

+ (x(3y 2 - x 2 ) - 2y 2 S y )c 4 w }D (3)}, (39) 

Ag r {x,y) = -AT{x~y), (40) 

Af{x,y) = -Af{x^y), (41) 
with B (i), Co(i), and D (i) the two-, three- and four-point scalar functions given in the Appendix. 
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III. NUMERICAL RESULTS AND DISCUSSION 



In the previous section the amplitudes for the sets of diagrams shown in Figs. ^]-^]have been arranged in an U c (l) 
gauge invariant form. We have shown that the respective amplitudes are SU(2) gauge independent and U c (l) gauge 
invariant for each set of diagrams in Figs. || and |3[ This property enables us to discuss separately each contribution. 
In particular, we would like to study the role played by the off-shell U (l) gauge structures of the Z boson and the 
neutrino. From (^) and ([n]) the double phase space integral for the decay width can be expressed as 

T(Z -> vvi) = [ dx [ dy{\MzA 2 + |AV| 2 + \M Box i\ 2 + \M Box5 \ 2 

Box4 + M Box5 Box 4 + M Bo.ro )Ml, 

+ M Boxi M Boxb ]}. (42) 

In order to disentangle the relative importance of the contributions arising from the ZZ*j and vv*^ couplings as 
well as the box diagrams, we have evaluated separetely the corresponding terms in (^|). In the case of the box dia- 
grams contributions, we have obtained separetely the contribution arising from each of the five U c (l) gauge structures 
defining these amplitudes. The result for the branching ratio B(Z — > vvy) will follow inmediately. 

A calculation of this type requires a careful numerical evaluation of the Passarino-Veltman scalar functions, in 
particular of those arising from diagrams with resonant effects, as it is the case for the graphs shown in Fig. |J. 
Since some of the involved scalar functions can not be expressed in terms of elementary functions, it is necessary to 
evaluate them by means of more elaborated methods, such as the Fortran program FF . We have made this task 
in different ways to get a reliable answer. In the case of the ZZ*^j and vv*^ couplings contributions, we get simple 
analytical solutions for the corresponding amplitudes, so we programmed the required integrals in Mathematica [23] ] 
and Fortran language together with the FF package. We have found an excellent agreement between the obtained 
results. In the case of the box diagrams, the four-point functions were analized by several methods. Explicit solutions 
of this functions were numerically evaluated and the results were compared with the corresponding evaluation using 
the FF package. In this way we are certain |M that we get the correct answer for these functions. We found delicate 
cancellations arising from the box amplitudes and in order to get the required accuracy different integration methods 
were used. 

Not all the box diagrams present troubles for their numerical evaluation. The most serious difficulties arise from 
the box diagrams shown in Fig. [|. In particular, the graphs (b) and (c) of this set induce collinear singularities. This 
fact implies, on one hand, that we must keep the lepton masses within the arguments of the scalar functions and, on 
the other hand, that we should have special care with the numerical evaluation of the corresponding amplitude. In 
order to get a reliable answer for these diagrams, we have used an entirely independent check. The idea is just to 
evaluate the amplitude via unitarity conditions 

ImM BoxA =^[d$ l M(Z^l + r)Mt(i?vj^l + r), (43) 

where is the phase space for the intermediate lepton pair. From the tree-level amplitudes of the right-hand side 
of this equation we can see that only the graphs |](b) and |(c) have collinear sigularities. The result for the left-hand 
side of (fig]), obtained by the above described method, agrees nicely with the one obtained for the right-hand side. 

We now proceed to discuss our results [J First, the coupling contribution to the decay width T[Z — > Dwy) is 

2.02 x 10~ 10 GeV. On the other hand, in the case of the ZZ*^ coupling only the top quark family gives an important 
contribution, namely 5.74 x 10~ 12 GeV, while the remaining contributions are negligible small. This fact can be 
understood because the ZZ*"/ coupling contributions have a strong dependence on the mass gap existing in each 
fermion family. In particular, an exact mass degeneration implies a zero contribution. This nondecoupling effect of 
chiral fermions is well known in the SM and it is due to isospin breaking. This constitutes an important effect in 
radiative corrections proportional to ml — m\ . In the same footing, it is interesting to study a possible enhancement 



3 The values used for the SM parameters were taken from the Particle Data Group |2f 
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for the ZZ*^f amplitude arising from a fourth family of chiral fermions with SM assignment of quantum numbers. 
The most promising scenario arises from a heavy top-like quark with a light bottom-like quark and a light tau-like 
lepton. In order to give an estimation of this contribution, we have used the LEP1 lower bounds for the tau-like 
lepton and the bottom-like quark: m' T < 80 GeV and m' b < 46 GeV; also, a variation of m' t in the 300-850 GeV 
range is assumed, where the upper limit for the top-like quark mass was taken from the Ap constraint over the fourth 
family mass difference |2(|. The corresponding contribution, as a function of Am = m' t — m' b shows an enhancement 
of almost 5 times with respect to the third family contribution, this enhancement is rapidly reached at Am = 300 
GeV. Finally, in the case of the box diagrams as well as the interference terms, their contribution is of the same order 
as that coming from the vv*^ coupling. We have found that the main contribution of the box diagrams arises from 
the U e (l) gauge structure corresponding to the neutrino MDT, whereas the U e (l) gauge structures given in ( |2^ ) lead 
to a contribution of order 10~ n -10~ 13 GeV. 

In summary, considering all the terms in ( |42| ) and assuming lepton universality we have B(Z — > Dvj) — 7.16 x 10~ 10 
for the three SM neutrino species. The main contribution comes from the neutrino MDT U c (l) gauge structure, while 
the Z boson EDT contribution is two orders of magnitude below. 



IV. FINAL REMARKS 

We have presented a complete calculation of the SM decay Z —> Dvj using a nonlinear gauge. We have found that 
this decay is characterized by the off-shell electromagnetic properties of the neutrino and the Z boson: a neutrino 
MDT and a Z boson EDT induced by the vv*^ and the ZZ*~f couplings, respectively. We found that the Z boson 
EDT contribution may be enhanced by a fourth generation of chiral fermions; however, this effect is not enough to 
increase significantly B(Z — ► vv~f). On the other hand, the main contribution to the decay Z — > vv^j arises from 
the neutrino MDT. This contribution depends strongly on the momentum transferred by the Z boson but it is not 
sensitive to a possible neutrino mass. Another interesting fact found in our calculation is that the box diagrams 
contributions can not be neglected because they induce a neutrino MDT U c (l) gauge structure with a form factor of 
the same order of magnitude as that coming from the vv**y coupling. It is important to notice that this is true even 
though there is no SU(2) nor U c (l) gauge link between the vv*^/ and box diagrams. This scenario was foreseen in a 
recent study of the decay Z — > vvf within the effective Lagrangian approach [pj: the operators inducing the ZZ*~f 
coupling are suppressed with respect to the dimension-six and dimension-eight operators inducing the vv**y coupling 
and the box diagram contributions. Therefore, it is expected that the dimension-six and dimension-eight operators 
contribute substantially to B(Z — > vvy), as it was the case in our calculation. Finally, we would like to point out that 
the decay Z — > Diry leads to a branching ratio of the same order of magnitude than the one expected for the radiative 
decay Z — > 777 [27] . In this case, the experimental bound for the respective branching ratio is also about four orders 
of magnitude below the SM prediction. Therefore, these decay modes of the Z boson should be of interest for people 
searching for deviations coming from new physics effects. 
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APPENDIX: PASSARINO— VELTMAN SCALAR FUNCTIONS 

In this appendix we present the arguments of the scalar functions involved in the Zvv^j interaction. We follow the 
notation of ||l9|| and the lepton mass was denoted by m;. 

Two-point scalar functions 

B (l) = B (xm 2 z ,m'i,rriw), 
.Bo (2) = B (ym%,mf,Tn% r ), 
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B (3) = B (m 2 z , m 2 w ,m 2 w ), 
B (4) = B (8m 2 z ,m 2 v ,m 2 v ), 
B (5) = B (m 2 z ,mf,mf), 

B (6) = B (6m 2 z ,mf,m 2 ). (1) 



Three-point scalar functions 



C (l) = C o (0, (5m|, m|, mfy, m^, m^), 
Cb (2) = C (0, 0, Sm%, m 2 w , mf, m 2 w ), 
C (3) = C (0, 0,xm|,TO^,m^,mf), 
Go (4) = Co(0,0,ym|,m^r,m^,rnf), 
Co (5) = C o (0,a;m|,m|,m^,mf,m^), 
C (6) = C o (0, ym|, m|, m^, mf, m^), 
C (7) = C o (0, Q,xm%,m 2 v ,m 2 ,m 2 ), 
C (8) = C o (0, 0, m 2 w ,mf, mf), 
C (9) = C o (0, xm%,m 2 z , mf,m 2 w , mf), 
C (10) = C (0,ym%,m% ) m 2 ,m 2 V) rni), 

C (ll)=C (0,0,6m 2 z ,mf,m 2 w ,mf). (2) 



Four-point scalar functions 

B (l) = B o (0, 0, 0, m|, xm%,Sm%, m\y, mfy, mf, m^), 
D (2) = D o (0, 0, 0, m|, ym|, <5m|, m^, m^, , m^r), 
B (3) = D (Q,0,Q,m 2 z ,xm 2 z ,ym 2 z ,m 2 v ,mf,mf,m 2 v ), 
B (4) = Du{Q,0,Q,m 2 z ,xm 2 z ,ym 2 z ,mf,m 2 v ,m 2 v ,mf), 
B (5) = D a {Q,Q,Q,m 2 z ,xm 2 z ,8m 2 z ,mf,mf,m 2 v ,mf), 

B (6) = B o (0, 0, 0, m|, m|, <5m|, mf, mf^, , ). (3) 

It is important to notice that lepton mass can be droped out from the above equations, except in Co (7), Bo (5) and 
Bo (6), where collinear singularities arise. 
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FIG. 1. Z* — 7 mixing term. Its contribution vanishes in the nonlinear gauge. 
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FIG. 2. ZZ*^f coupling contribution. The crossed diagram must be added. 
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